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Q ■ Abstract 

^vq I We have searched for hadronic axions which may be produced in the Sun by a 
bremsstrahlung-hke process, and observed in the HPGe detector by an axioelectric 

^. effect. A conservative upper Hmit on the hadronic axion mass of rua ^ 334 eV at 

D I 95% C.L. is obtained. Our experimental approach is based on the axion-electron 

pi i' couphng and it does not include the axion-nucleon coupling, which suffers from the 

^ . large uncertainties related to the estimation of the flavor-singlet axial- vector matrix 

i__i| element. 
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Q ■ The axion, a light pseudoscalar particle associated with the spontaneous break- 
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ing of a f/(l) Peccei-Quinn symmetry [T], was introduced to explain the ab- 
sence of CP violation in strong interactions. The mass of the axion satisfies 
rUa = 6eV x 10^ GeV/ fa [2]- The symmetry-breaking scale (or the axion de- 
r> I cay constant) fa, however, is left undetermined in the theory. The present 

c^ ■ astrophysical and cosmological considerations [2] have placed bounds on the 

parameter, and are consistent with 10~^ eV < rria < 10~^ eV. At the lower 
end of this constraint, axions are a viable cold dark matter candidate, and ex- 
perimental attempts to detect their presence are in progress [3]. Besides this 



range of allowed axion masses, for hadronic axiono [1] there exists another 
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^ Although hadronic axions do not couple directly to ordinary quarks and charged 
leptons, their coupling to nucleons and electrons is not zero due to the axion-pion 
mixing and radiatively induced coupling to electrons. Their coupling to photons is 
model dependent. It was shown by Kaplan [5] that is possible to construct hadronic 
axion models in which the axion-photon coupling is significantly reduced and may 
actually vanish. 
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window of 10 eV < nia < 20 eV [2], as long as the axion-photon coupling is suf- 
ficiently small [5]. These bounds arising from arguments concerning the super- 
nova 1987A cooling and axion burst, however, are model dependent and with 
large uncertainties. In respect of the early universe, axions in this hadronic 
axion window can reach thermal equilibrium before the QCD phase transition 
and hence, like neutrinos, they are also candidates for hot dark matter [6ff7] . 



As axions couple to photons, electrons and nucleons, the Sun would be a strong 
axion emitter. Hadronic axions of continuous energy spectrum, with an aver- 
age energy of 4.2 keV, could be produced abundantly in the solar core by the 
Primakoff conversion of thermal photons in the electric fields of charged parti- 
cles in the plasma. The ongoing CAST experiment at CERN searches directly 
for these axions by pointing a decommissioned Large Hadron Collider proto- 
type magnet (with a field of 9 T and length of 9.26 m) toward the Sun. Thanks 
to the powerful magnet and the installed x-ray focusing mirrors (which reduce 
noise), the detection sensitivity is considerably improved relative to previous 
experiments of this kind [8]. The obtained upper limit on the axion-photon 
coupling of 8.8 x 10~^^ GeV~^ [5] also supersedes, for a broad range of ax- 
ion masses, the previous limit derived from energy-loss arguments on globular 
cluster stars [2] . In the case of hadronic axions with strongly suppressed photon 
couplings the Primakoff rate is negligible. Some nuclear processes have been 
proposed as sources of solar monoenergetic axions [TOl[TTj and experiments 
based on the detection of hadronic axions with suppressed axion-photon cou- 
plings were reported by several authors |12fl3lll4lll5fl6fl7lll8fl9] . To date the 
best experimental limit on mass of these axions is set to be around 216 eV 
p!7] . However, as pointed by the author himself, this experimental approach 
based on the axion-nucleon coupling suffers from a poorly constrained flavor- 
singlet axial-vector matrix element that affects the axion-nucleon interaction 
strongly. For example, if a recent value of this matrix element of ^ 0.3 [2U] is 
taken into account, the obtained result becomes weaker, rria < 515 eV (95% 

c.L.) m- 



Here we focus our attention to emission of the hadronic axions via their ra- 
diatively induced coupling to electrons [2T] in the hot solar plasma and only 
the bremsstrahlung-like process, e~ + Ze{e~) — > e" + Ze{e~) + a, is used 
as a source of axions. In addition, following the calculations in [22], we have 
estimated that the contributions of electron-electron collisions to the axion 
bremsstrahlung emission from the Sun are negligible in the 2.0 - 3.8 keV en- 
ergy region which is of interest in our experiment. Therefore, only scattering 
of electrons on protons and He nuclei was considered. Using theoretical pre- 
dictions of Zhitnitsky and Skovpen [23] for the axion bremsstrahlung due to 
electron-nucleus collisions, for the case where the Born approximation is valid 
and E ^ m^ {E is the total energy of axion), we can write the differential 



solar axion flux at the Earth as 

Rq 



1 ^ 



dE ~ iird'c:. 

^ 



Here ^0 is the average distance between the Sun and the Earth, Rq denotes 
the solar radius while A'^h and iVne are the number density of hydrogen and 
helium nuclei in a given spherical shell in the solar interior at the radius r, 
respectively. The Maxwellian distribution at the temperature T{r) for the 
nondegenerate, nonrelativistic incident electrons of velocities Ve and kinetic 
energies Tg is given by 






^e= v., ^,3/2 ^e, (2) 



where A'^e is the number density of the electrons at the radius r, and k is the 
Boltzmann constant. The differential cross section for the axion bremsstrahlung 
process (for Z=l) [23] is designated by daa/dE. Integrating the expression of 
Eq. ([1]) over the BS05 standard solar model [23] we flnd that the expected 
solar axion flux at the Earth is $„ = g^^ 1.3 x 10^^ cm~^ s^^ with an average 
energy of 1.6 keV and an approximate spectrum 

^ = gl 1.55 X 10^5 (Ekev)°-^^ e-^-v cm'^ s'' keV"^ , (3) 

dE 

as shown in Fig. [H Here -Ekcv = E/keV. The corresponding solar axion lumi- 
nosity is calculated to be La = gle'^A x lO^^L©, where L© = 3.84 x 10^*^ W 
is the solar photon luminosity. In particular, the coupling exploited in our 
experiment, of the hadronic axion to electrons, is given numerically by [1] 



gae =6.6 X 10 ^^mev 



E 

— (23.2-lnmev) -14.8 



(4) 



where E/N is the model-dependent ratio of electromagnetic and color anoma- 
lies, and rricv = rna/eV. For hadronic axions that have greatly suppressed 
photon couplings {E/N ^ 2) j5] this is gae ~ (2.0 -^ 1.2) x 10"^^mev over 
a broad range of axion masses of 1 ^ 1000 eV. We note that astrophysical 
arguments related to stellar evolution {La < Lq) [23] and helioseismology 
{La < 0.2 Lq) [26] would imply upper bounds on the axion-electron coupling 
of 6.5 X 10"^^ and 2.9 x 10^^^, respectively. By using Eq. (j4]) they translate to 
upper bounds on the hadronic axion mass of 515 eV and 210 eV. 

In this Letter, we report results of our search for solar hadronic axions which 
could be produced in the Sun by the bremsstrahlung-like process and detected 
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Fig. 1. Differential solar axion flux at the Earth, derived by integrating Eq. ([1]) 
over SSM [24] up to r = Rq (red line), r = 0.2 Rq (blue line), and from r = 0.2 Rq 
to r = Rq (light blue line). The axion-electron coupling Qae is defined in Eq. ^. 

in the single spectrum of an HPGe detector as the result of the axioelectric 
effect on germanium atoma ^ I. The x rays accompanying the axioelectric ef- 
fect will be subsequently absorbed in the same crystal, and the energy of the 
particular outgoing signal equals the total energy of the incoming axion. Be- 
cause in this experimental set-up the target and detector are the same, the 
efficiency s of the system is substantially increased, i.e., e ^ 1. The HPGe de- 
tector with an active target mass of 1.5 kg was placed at ground level, inside 
a low-radioactivity iron box with a wall thickness ranging from 16 to 23 cm. 
The box was lined outside with 1 cm thick lead. The crystal was installed in 
a standard PopTop detector capsule (Ortec, model CFG/PH4) with a beryl- 
lium window of 0.5 mm. The HPGe preamplifier signals were distributed to 
a spectroscopy amplifier at the 6 /is shaping time. A low threshold on the 
output provided the online trigger, ensuring that all the events down to the 
electronic noise were recorded. The linearity and energy resolution have been 
studied by using various calibrated sources and, in particular, in the lowest- 
energy region mainly a ^^^Am source (13.9 keV x-rays and their escape peak of 
3.9 keV). Data were accumulated in a 1024-channel analyzer, with an energy 
dispersion of 63.4 eV/channel, in 20-hour cycles with total time of collection of 
275 days. In long-term running conditions, the knowledge of the energy scale 
is assured by continuously monitoring the positions and resolution of the In 
x-ray peaks which are present in the measured energy distribution. Drifts were 
< ±1 channel and a statistical accuracy of better than 0.6% per channel was 
attained. Thus, all of the spectra were summed together without applying any 
gain shifts or offsets and are shown in Fig. [2l Energy resolution (FWHM) was 
estimated to have a value of 660 eV at the photon energy in the 2.0 - 3.8 keV 
region, where the axion signal is expected. The former bound is imposed as 



^ Using theoretical predictions of [23] we estimate that for ~keV axions the axio- 
electric effect is about three orders of magnitude stronger than the axion-to-photon 
Compton process. 
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Fig. 2. Total measured energy spectrum showing also x-ray peaks from various 
materials. 

the analysis threshold due to the electronic noise while the latter one is due 
to the energy distribution of the solar axions. 

Our experiment involves searching for the particular energy spectrum in the 
measured data, 



-C'n + l OO ,^ 

dE' j -^aae-.eR{E\E)dE 



(5) 



produced if the solar axions are detected via axioelectric effect. Here 5"^ is 
the number of counts detected in detector energy channel n, 7VGe=l-24xlO^^ 
is the number of germanium atoms in our 67 mm in diameter and 80 mm 
thick HPGe crystal, and t is the time of measurement. Factor 2 is number of 
electrons in (Li-) Mi-shell while their binding energies are -E'b(Li)=1.413 keV 
and i?b(Mi)=0.181 keV. The axion response function of the detector R{E , E) 
is well represented by a Gaussian (to describe full energy peak) with the 
amplitude, normalized on the efficiency e, and width of 0.42 FWHM. Using 
Eq. (3) from [13], the cross section for axioelectric effect per electron was 
calculated including contributions from 2s and 3s initial state electrons and 
may be written in the ma -^ Sz E <^ rrif, approximation and for E > Ei,(Li) 
as: 



(Jr, 
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Here crxh = 6.65 x 10~^^ cm^ is the cross section for Thomson scattering, 
a=l/137 is the fine-structure constant, Z=32 is the atomic number of Ge, 
and rrie is the electron mass. Equation ([6]) is consistent with a more general 
expression aae^e = (oae/a) (-E/2me)^ cTph [27], where aae = qL/^'^ and cTph is 



the photoelectric cross section. 

We have used the fit method to analyse the data in order to find upper limits 
on Qae {fna)-, cis described below. Similar approaches have been used by other 
groups [28] to extract (conservative) upper limits in a case like this where 
direct background measurement is not possible (the Sun cannot be switched 
off) and the signal shape is a broad smooth spectrum on top of an unknown 
background spectrum. 

In the fit method we have used the following procedure to make the estimation 
of an axion signal. The experimental data in the energy interval 2.0 - 3.8 keV 
(corresponded to MCA channels of 32 to 60) was fitted by the sum of three 
functions: the electronic noise expectation (T„), the background expectation 
{Bn) and the effect being searched for {Sn)- To estimate the electronic noise 
the data were fit to a four-parameter function in the region below 2 keV 
(corresponded to the channels of 4 to 24); it has the form T„ = An^ex^iCn^), 
where A = 1.46 x 10^ B = 6.85, C = -0.858, and D = 1.08, see Fig. [3] (top 
panel). As a simple background estimation, suitable for the present purposes, a 
first-order polynomial Bn = a + bn has been assumed. Values of a = 4.54 x 10^ 
and b = —42.5 were found by fitting Bn with the data in the region above 3.8 
keV, corresponded to the channels of 70 to 110. From the current limits on 
the axion mass [T7fT8] we estimate that the axion signal is negligible in this 
channel interval. Contribution of axions was then investigated by extrapolating 
the electronic noise T„ and background B^ in the 2.0 - 3.8 keV region, where 
the detection of axion events should be the most efficient. Inserting Eqs. ([3]) 
and Qj in Eq. ([5]), one can express the expected axion spectrum Sn as a 
function of gj^^. As the best values of A, B, C, D, a and b had been already 
determined in the regions where axion contributions could be neglected, only 
(jf^g was varied in the x^ comparison. We were using g^^ instead of gae as the 
minimization parameter because the signal strength (i.e., number of counts) 
is proportional to g^^. The results of the analysis are consistent with g^^ = 
(3.6 ± 0.1) X 10"^^ at lo" level. Figure [3] (bottom panel) displays the results of 
our fit. Since we cannot exclude more complicated scenarios in which both the 
background and the electronic noise modify significantly the energy spectrum 
of signal events in the region of 32 to 60 channels, we treat the residuals 
from the noise and background expectation conservatively as an upper limit 
on the axion signal. In other words, we adopted a criteria commonly used in 
experiments where direct background measurement is not possible |28] that 
the theoretically expected signal cannot be larger than that observed at a given 
confidence level. In this way we obtain the standard 95% C.L. upper limit [2] 
Qae ^ 4.4 X 10^^^ which translates through Eq. (jlj) to m^ < 334 eV. Because 
the precise form of the background expectation is not known, we tried also 
other plausible forms, such as Bn = a + bn + cn^ and Bn = a + bn + cn~^. The 
former yields g^^ = (2.5 ± 0.2) x 10"^^ ^suiting in m„ < 310 eV (95% C.L.) 
while the latter provides g'^^ = (1.4 ± 1.9) x 10~^^ implying the more stringent 
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Fig. 3. Top panel: low-energy data are shown together with the best fit for the 
electronic noise expectation Tn (solid line) and its extrapolations (dashed line). 
Bottom panel: residuals from the noise and background expectation are shown to- 
gether with the expectations for best fit gae = 4.4 x 10~^^ (red line) as well as for 
gae = 5.0 X 10~^^ (blue line) and gae = 3.7 x 10~^^ (green line). On both panels 
arrows indicate the energy region of interest for axion signal. 

limit rria ^ 190 eV (95% C.L.). One can see that subtraction of the background 
expectation in a form oi a + bn + cn~^ cancels the excess of events and the 
obtained result is very close to the sensitivity of the experiment. Namely, if 
we assume that the measured spectrum in the region of interest for axion 
signal is compatible with the background expectation, the fluctuation due to 
statistical uncertainty imposes a limit on the maximum allowable number of 
axion events. In this case (the mean is zero and statistical uncertainty for g'^^ 
is cr = 2 X 10~^^) we obtain Qae ^ 2.4 x 10~^^ corresponding to nia < 170 eV 
at the 95% confidence level. 



In conclusion, our measurements based on the coupling of the hadronic axions 
to electrons set a conservative upper limit on the axion mass of rua ^ 334 eV 
at 95% confidence level. This upper limit is comparable to the current limits 
[T71IT8] which are extracted from the data taken in the experiments based on 



the axion-nucleon coupling. We point out the derived hmit is free from the 
large uncertainties and ambiguity, associated with the estimation of the flavor- 
singlet axial- vector matrix element, which are related to experiments based on 
the axion-nucleon coupling. New experiments with lower electronic noise and 
reduced background could improve our search for the hadronic axions. Note 
that the axioelectric absorption peak in germanium occurs at lower incoming 
axion momentum when E coincides with the binding energy of atomic shells. 
At these energies, atomic bound state effects lead to large enhancements in 
the detection rates of axions, similarly to those in the photoelectric effect. For 
example, the expected number of axion events for E > lA keV is about 3 
times that for E' > 2 keV. Thus an improvement in the discovery potential 
may come from examination of the Ge data for energies near the Li-shell 
peak of 1.413 keV. Upgradings of the set-up to reduce the electronic noise 
and permit lowering of the energy threshold below 1 keV are already foreseen 
and in preparation. In the near future we expect that measurements with the 
array of 10 - 20 ultra-low-energy IGLET germanium detectors with an anti- 
Compton shield could improve significantly our understanding on hadronic 
axions. 
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